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Method for Determining Positions of Targets by Bistatic Measurements Using 
- Signals Scattered by the Targets . 

The present invention relates to a method for determining positions of targets by 
bistatic measurements using signals scattered by the targets. Also the velocities of 
the targets can be determined. The method comprises a rapid bistatic association 
method which is suitable for, for instance, a network of radar stations in the manner 
of AASR (Associative Aperture Synthesis Radar) although there may be further 
fields of application. AASR is described, inter alia, in Swedish Patent 0101661-7, 
which is herewith incorporated by reference. In the following, the description will be 
concentrated on the new method of associating by bistatic measurements only. 

First the fundamental problem that is solved by the invention will be presented. N s 
stations (for instance radar stations) are imagined to be set out in the space (R 3 ). 
The stations are designated Sj, j=1, ...N s and their position vectors are designated 
p 3 j=1, ... N s . In addition to the stations, there are also N t moving targets which are 

to be detected. They are designated tj, i=1, ... N t and have corresponding time- 
dependent position vectors n = r^t), i=1, ...N t . 

Each station is capable of measuring distances (up to a certain maximum distance) 
and radial speed for each target. Thus, the station Sj, 1 < j < N s will, at a certain point 
of time, measure 

dj(k) = |r k - Pj |, k=1,2, ... N dj <N { 

Vj (k) = (d/dt)|r k -pj |, k=1, 2, ... N dj < N t 

For stations that are sufficiently close to each other, also bistatic measurement infor- 
mation is obtained, i.e. transmitting from one station and registration at another. For 
the pair of stations (Sj, Sj) it means that the following is registered 

dy(k) = |r k - Pi |+|r k - Pj |= d { (k)+ dj(k), k=1 , 2, . . . N dlj < N t 

Vij(k) = (d/dt)|r k - Pj |+(d/dt)|r k - Pj |= v { (k)+ Vj(k), k=1 , 2, N dij < N t 



It is to be noted that with these designations, d n (k)=2di(k), v H (k)=2v j (k) > i=1,2 

k=1,2, ... 
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For each sensor (monostatic or bistatic geometry) targets are thus registered in 
respect of distance and Doppler. It is a priori not possible to know which registration 
from one sensor is associated with a certain registration from another sensor, i.e. 
originating from the same target. If registrations from different sensors are paired 
5 incorrectly, false targets, ghost targets, arise. The problem of association is to 
discriminate, among all conceivable possibilities of combining sensor data, 
corresponding to conceivable target candidates, between correct combinations 
(targets) and false combinations (ghosts). 

10 The maybe most straight-forward method is to consider three neighbouring stations 
and their monostatic registrations, which for the sake of simplicity are assumed to be 
N in number. These measurements can be combined in N 3 ways where each combi- 
nation corresponds to a target position which is determined up to reflection in the 
plane containing the three stations. (Certain combinations can be incompatible, 
15 corresponding to false candidates.) These ~N 3 candidates can then one by one be 
■*v: compared with the bistatic measurements and either be rejected or accepted. The 
s problem of such a method is that it will be very slow if the number of targets, N, is 
* large. For this reason, more efficient association algorithms have been developed. 

20 Each target is to be determined in respect of position as well as velocity, i.e. they are 
to be positioned in a six-dimensional state space. The number of cells in the state 
space can be very large (~ 10 18 ), which means that traditional projection methods 
will be irretrievably slow. 

25 The above Swedish Patent 0101661-7 discloses a method of attacking the problem 
of association by designing a sensor network so that each target is registered by 
many sensors (monostatic and bistatic), i.e. a high degree of redundancy is obtained 
in the system. Then the state space is divided into a manageable number of 
relatively large cells. 

30 

If the cells are just large enough, it will be possible to reject many of them, i.e. they 
cannot contain any targets, for the following reasons. If the cell contains a target, all 
(or almost all) of the possible sensors that can register targets in the current cell, will 
indicate such a registration. On the other hand, if the cell is empty, some sensors, 
35 .and yet not too many, will still indicate registrations (from other targets) that are 
compatible with the cell in question. Owing to redundancy, a sufficient number of 
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sensors will indicate the cell as empty, and it can be rejected. When the number of 
cells is thus reduced, the surviving cells are divided into smaller cells and the 
process is repeated. The process is repeated until the cells in the state space have 
reached the desired size. As the cells are becoming smaller, fewer and fewer ghost 
5 targets will survive, so that, when interrupting the process, practically only real 
targets are left. What speaks in favour of this method is that it uses (but also 
requires) the redundancy of the sensor network. However, it is not yet quite clear 
how rapid the method may eventually be. 

10 An alternative method is disclosed in Swedish Patent 0101662-5, which is herewith 
incorporated by reference, and implies that use is made of certain symmetries of the 
' combination sensors - measurement data. Given two stations, the two monostatic 
measurements, together with the bistatic measurement, will share a symmetry, viz. 
that the three measuring geometries are all insensitive to rotation of the targets 
15 about the axis extending through the two stations. This means that it is possible to 
remake an initial rapid screening of the candidates and delete a large number of false 
^associations (ghosts). The subsequent final association will then be significantly 
i. more rapid. A drawback, however, is that the monostatic measurements will be 
important, which may be disadvantageous in connection with reconnaissance of 
20 Istealth targets. 

The present method according to the invention is based on using a rapid method 
Where only bistatic measurements are utilised. Furthermore the method manages a 
certain dropout of sensors in a better way than the method that has been discussed 
25 directly above. The method solves the current problem of association by being 
designed in the manner as is evident from the independent claim. Advantageous 
embodiments of the invention are defined in the remaining claims. 

Before a more detailed description of the invention, first a multistats network of 
30 ground radar stations will be contemplated, in which each radar station transmits 
radar pulses that are scattered towards flying targets and are then received by the 
surrounding stations. There will then be a situation involving a large number of 
bistatic measurements (i.e. the transmitting and the receiving station are located in 
different positions) and also monostatic measurements which, however, are not 
35 used in the invention. The bistatic measurements contain information about the total 
distance transmitter-target-receiver and corresponding Doppler information. A 
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coherent air situation image is then to be created from all these measurements. This 
problem, the problem of association, is non-trivial if there are a large number of 
targets. 

5 For intuitive understanding of the invention, a simple case is taken into considera- 
tion, involving only one target, m 1t and four stations s 1f s 2 , s 3l s 4 . Now assume that 
the measurements d 12 , d;*, d 13( d 24 are performed, where d 0 means the total 
distance Sj - nr^ - Sj. It will be appreciated that d 12 + d 34 = d 13 + d 24 must be the case 
since both expressions mean the total distance from the target to the four stations. 
10 . 

If there are now N targets instead, the above observation can be used to correctly 
associate data in the following manner. All conceivable combinations of data of the 
type d 12 and d 34 are formed; they will be N 2 in number. In the same way, N 2 combi- 
nations of data of the type d 13 and d 24 are formed. These combinations are sorted 
15 and compared, and only. sums from the two amounts that are equal (within a given 
tolerance) can correspond to real targets. The same discussion can be used about, 
the Doppler velocities which thus give a further screening. In this way, quick and 
easy association of measurement data:can be effected. 

In general, the transmitters and receivers must be positioned and the range of the 
transmitters must be chosen so that a target at an arbitrary point within the position 
space can be measured via scattering in the target of at least four cooperating 
: bistatic pairs of transmitters and receivers. The number of transmitters and receivers 
can be large. At least four such cooperating pairs are selected among these bistatic 
pairs- to perform the association and the determination of the distance. 

Below follows a more systematic presentation of the calculations: In order to obtain a 
simple description, the following (non-critical) assumptions are made. Assume that 
there are four stations and N targets, which all are seen by all sensors (monostatic 
as well as bistatic). 

Input data is thus (monostatic measurements) 



dj(k) = |r k - Pj |, k=1, 2, ... N , j=l f 2,3,4 

Vj(k) = (d/dt)|r k - Pj |, k=1,2, ... N, j=1,2,3,4 
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and (bistatic measurements) 



dij(k) = |r k - Pi |+|r k - Pj |= di(k)+ dj(k) 



k=1, 2, ... N, j=1, 2,3,4 
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Vij(k) = (d/dt)|r k - pj |+(d/dt)|r k - Pj |= v s (k)+ Vj(k), k=1,2, ... N, j=1 ,2,3,4 

where i = j in the bistatic case corresponds to monostatic measurements, i.e. i * j 
can be assumed if desirable. 

It is to be noted that, for instance, for station j, with the monostatic measurement 
dj(k), k=1,2, ... N, it is not possible to know which measurement belongs to a certain ' 
target, i.e. the measurements are to be regarded as a set which is as a suggestion 
sorted according to distance. In this way, there is no connection between a certain 
index k which belongs to two different sensor registrations. 

The method is now based on the following observation: For each registered target 
- (not candidate, but real target) there must be a k, a k\ an I and an l\ all between 1 
and N so that 

di2(k)+d 3 4(I)=d 13 (k , )+d 2 4(r) 
For the same k, k\ I, I', the following is also applicable 
Vi2(k)+v 34 (l)=v 13 (k>v 24 (l') 

The reason is that if the target has the space vector r t , it is applicable for the target 
that 

di2(k)+d 3 4(!)= |rt-PiMrt-p 2 |+|rt^p 3 |+|r t -p 4 | 
and in the same way that 

d 13 (k)+d 24 (])= |r t -pi|+|rt-p3|+|rt-p 2 |+|rt«p 4 | 



35 



so that they are equal. The argument for the velocities is identical. The suggested 
method now is as follows. 
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Step 1. Form the N 2 sums 

di2(k)+d34(l), 1 < I, k < N 
5 Sort them according to the total distance and designate them 
di 2+ 34(m), 1 < m < N 2 
Step 2. Proceed in the same way with 

10 

di3(k , )+d 24 (l'), 1 < I", k' < N 
so that the following will also be obtained (sorted) 
15 d, 3 + 24 (m'), 1 < m' < N 2 

Step 3. Associate targets from { d 12+3 4(m)} m=1 , 2 ... N 2 with targets from 

■ { di 3 +24(m , )}m4 S i^„. N 2 if 

a: 

20 | d 12+3 4(m)- d 13+2 4(m')| < suitable tolerance 

Step 4. Investigate, and keep associated targets if they also satisfy 
I v 12+3 4(m)- v 13+24 (m')| < suitable tolerance 

25 

"Suitable tolerance" in Step 3 is determined, inter alia, by the transmitted signal 
bandwidth, the purpose of the processing and hypotheses about size and number of 
the targets. Usually it is from about one meter to some twenty or thirty meters. 
Correspondingly, "suitable tolerance" in Step 4 is usually a few meters/second. 

30 

To establish that this really results in a rapid method, the following rough estimate 
may be used. Assume that there are many targets so that they are of the same 
magnitude as the number of distance bins and the number of Doppler bins. This 
common number is again designated N. It may then be estimated that, since the total 
35 number of cells (= number of distance bins by the number of Doppler bins) is the same 
as the number of candidates in for instance { d 12+34 (m)} m=1 , 2 ... N 2 , each such candidate 
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will be paired with typically a false candidate from { d 13 +24(m , )} m - 1 , 2 ... n2. The number of 
candidates according to the above procedure thus is ~ N 2 (fewer with fewer targets), 
which is a great reduction compared with N 3 . Further processing can then take place 
by comparing with the remaining bistatic geometry {d 14+2 3(m ,, )}m'tei l 2... N 2 , the mono 
5 static measurements or measurements involving other stations. 

It is also to be noted that it is possible to involve {d 14+ 23(m")} m » s1i 2... N 2 from the 
beginning. This gives a possibility of having a redundancy, i.e. a possibility of 
managing a certain dropout in registrations, in the following way. The condition that | 

10 d 12+34 (m)- d 13+24 (m')| < "suitable tolerance" can be seen as if both d 12+34 (m) and 
d 13 +24(m') are to be close to a certain given value. By requiring instead that two of 
di2+34(m) l d 13+ 24(m') and d 14+23 (m") should be close to the indicated value (for some 
values of m, m' and m") there will still be a discrimination between false candidates 
(ghosts) and targets. However, it may be accepted that one of the measurements 

15 drops out. 

The calculations in their entirety require 0(N 2 log N) operations, and there are ~ . 
simple methods of really obtaining position and velocity from the candidates, i.e. 
after processing four bistatic distances are known for a certain candidate as follows: 

20 

Ir-p^ + |r-p 2 | = d 12 

|r - Pal + |r - p 4 | = d 34 
|r-Pi| + |r-p 3 | ■= dia 

|r - p 2 | + |r - p 4 | = d 24 

25 

It is, of course, interesting to know the value of r (position of the target), i.e. a 
method of solving the above system of equations. (p j5 i=1, 2,3,4, are the known 
positions/ position vectors of the stations and d 12 , d 34 , d 13 , d 24 are the measured 
r bistatic distances.) Generally seen, intersections of ellipsoids cause relatively com- 
30 plicated algebraic systems of equations, but in this case the system of equations can 
be solved by simpler methods. 

If the system of equations is regarded as a 4x4 system, it is obvious that it is dege- 
nerated. At the same time the condition d 12 -+ d 34 = d 13 + d 24 guarantees that there is 
35 a parameter solution. By selecting the origin of coordinates in p 4 so that 
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|r- p 4 | = |r| = r and introducing r as a parameter, the following equations are 
obtained 

I r-p! | = d 12 - d 24 + r 
5 | r-p 2 | = d 24 - r 

I r-p 3 1 = d 34 - r 

It is here possible to square the three equations, in which case r 2 can be deleted, 
and obtain the following (for some a.b.c.cc.p.y) 

10 

r • pi = ar + a 
r • p 2 = br + (3 
r • pi = cr + y 

15 The latter system of equations can then be solved in a fairly straight-forward way. 
However, there will be two different cases in dependence on whether {pi}j=i )2(3 is 
linearly dependent or not. 

The case of the velocities, is similar, the following system of equations will be 
20 obtained 

u^-v + u 2 -v 

u 3 V + U A V 
• v + Q 3 - V 

25 u 2 -v + u A -V 

where u } = p: — z-j , v = r , / = 1 ,2,3,4. The system of equations can be processed 
T in the same fundamental way as the previous system of equations/ 

30 The invention can be implemented in high-level languages which are suitable for 
calculations, such as MatLab, C, Pascal, Fortran etc. 



= ^12 
= ^34 

= ^24 



WO 03/107033 




PCT/SE03/00988 



Claims: 

1. . A method for determining positions of targets in a position space using 
signals scattered by the targets, comprising use of a number, spread in known 
points in the position space, of transmitters and receivers of electromagnetic or 
acoustic signals, each bistatic pair of transmitter and receiver being referred to as a 
measuring facility, further comprising analysis of received signals, which includes 
determining of moments of transmission and reception according to generally 
accepted principles of radar and parameterisation of received signals as a function 
of the path of propagation between transmission point and reception point, but 
^without the conventional requirement in radar for directional information, the 
positions, being primarily determined by selecting the position of the transmitters and 
-. receivers and the range of the transmitters so that a target at an arbitrary point 
within the position space can be measured by scattering in the target by at least four 
•cooperating measuring facilities; by selecting an even number of cooperating 
measuring facilities, however at least 4, for the determination; by associating targets 
by calculating, in two independent ways, two sets of sums of distances between 
transmission points and targets and, respectively; targets and reception points, 
'based on bistatic distances, measured via the targets, for selected measuring 
facilities, sorting said two sums with respect to the distance, comparing these with 
.each other and establishing that the sums, calculated in the two different ways, 
which correspond with each other, while considering a margin of error that has been 
determined in advance, are stated to correspond to conceivable targets; and by 
calculating the positions of the targets from a system of equations for the bistatically 
measured distances, characterised by improving and completing the 
association of targets by performing calculations for bistatically measured Doppler 
velocities, corresponding to calculations for distances, and establishing that the 
Sums, calculated in the two different ways, which correspond with each other, while 
considering a margin of error that has been determined in advance, are stated to 
correspond to targets. 

2. A method as claimed in claim ^characterised by improving 
the association of targets by calculating the sum of all distances between the targets 
and the transmission points and the reception points, respectively, in a third way as 
the sum of bistatically measured distances via the target for measuring facilities, 
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sorting the sum with respect to the distance, comparing this with previously calcu- 
lated sums of distances, and establishing that the cases where two of the three 
sums, calculated in said different ways, correspond with each other, while consider- 
ing a margin of error that has been determined in advance, are stated to correspond 
to targets. 



3. A method as claimed in claim 1, characterised by improving 
the association of targets by calculating the sum of all Doppler velocities between 
the targets and the transmission points and the reception points, respectively, in a 
third way as the sum of bistaticaliy measured Doppler velocities via the target for 
measuring facilities, sorting the sum with respect to Doppler velocity, comparing this 
with previously calculated sums of Doppler velocities and establishing that the cases 
where two of the three sums, calculated in said different ways, correspond with each 
other, while considering a margin of error that has been determined in advance, are 
stated to correspond to targets. 

4. A method as claimed in claim 2 or 3, characterised by : 
improving the association of targets by requiring that all three sums, calculated in 
said different ways, correspond with each other, while considering a margin of error 
that has been determined in advance, for targets to be indicated. 

5. A method as claimed in claim 4, c h a r a c t e r i s e d by calculating 
the velocities of the targets from a system of equations for the bistaticaliy measured 
Doppler velocities. 

6. A system for determining positions of targets in a position space using 
signals scattered from the targets, comprising a number, spread in known points in 
the position space, of transmitters and receivers of electromagnetic or acoustic 1 
.signals, each bistatic pair of transmitter and receiver being referred to as a measur- 
ing facility, further comprising analysis equipment for storing and analysing received 
signals, which includes determining of moments of transmission and reception 
according to generally accepted principles of radar and parameterisation of received 
signals as a function of the path of propagation between transmission point and 
reception point, but without the conventional requirement in radar for directional 
information, the positions being primarily determined by the position of the trans- 
mitters and receivers and the range of the transmitters being selected so that a 
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target at an arbitrary point within the position space can be measured by scattering 
in the target by at least four cooperating measuring facilities; by the analysis equip- 
ment selecting an even number of cooperating measuring facilities, however at least 
4, for the determination; by the analysis equipment associating targets by calculat- 
5 ing, in two independent ways, two sets of sums of distances, between transmission 
points and targets and, respectively, targets and reception points, based on bistatic 
distances, measured via the targets, for selected measuring facilities, sorting said 
two sums with respect to the distance, comparing these with each other and 
establishing that the sums, calculated in the two different ways, which correspond 

10 with each other, while considering a margin of error that has been determined in 
advance, are stated to correspond to conceivable targets, and by the analysis 
equipment calculating the positions of the targets from a system of equations for the 
bistatically measured distances, characterised in that the analysis 
equipment improves and completes the association of targets by performing caicu- 

1 5 lations for bistatically measured Doppler velocities, corresponding to the calculations 
for the distances, and establishing that the sums, calculated in the two different 
ways, which correspond with each other, while considering a margin of error that has 
been determined in advance, are stated to correspond to targets. 

*- 

20 7. A system as claimed in claim 6, characterised in that the " 
analysis equipment improves the association of targets by calculating the sum of all 
distances between the targets and the transmission points and the reception points, 
respectively, in a third way as the sum of bistatically measured distances via the 
target for measuring facilities, sorting the sum with respect to the distance, compar- 

25 ing this with previously calculated sums of distances, and establishing that the cases 
where two of the three sums, calculated in said different ways, correspond with each 
other, while considering a margin of error that has been determined in advance, are 
stated to correspond to targets. 

30 8. A system as claimed in claim 6, characterised in that the 

analysis equipment improves the association of targets by calculating the sum of all 
Doppler velocities between the targets and the transmission points and the reception 
points, respectively, in a third way as the sum of bistatically measured Doppler 
velocities via the target for measuring facilities, sorting the sum with respect to 

35 Doppler velocity, comparing this with previously calculated sums of Doppler 

velocities and establishing that the cases where two of the three sums, calculated in 
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said different ways, correspond with each other, while considering a margin of error 
that has been determined in advance, are stated to correspond to targets. 

9. A system as claimed in claim 7 or 8, characterise d in that the 
analysis equipment improves the association of targets by requiring that all three 
sums, calculated in said different ways, correspond with each other, while consider- 
ing a margin of error that has been determined in advance, for targets to be 
indicated. 

10. A system as claimed in claim 9, characterised in that the 
analysis equipment calculates the velocities of the targets from a system of 
equations for the bistatically measured Doppler velocities. 




Abstract: 

The present invention relates to a method for determining the positions of targets by 
bistatic measurements using signals scattered by the targets. Also the velocities of 
the targets can be determined. The range of the transmitters is selected so that a 
target at an arbitrary point can be measured by scattering in the target by at least 
four cooperating measuring facilities. 

First the targets are associated. This occurs by calculating, in two independent 
ways, two sets of sums of distances between transmission points and targets and, 
respectively, targets and reception points. Subsequently, said two sums are sorted 
with respect to distance, compared with each other, and the sums that correspond 
with each other, while considering a margin of error that has been determined in 
advance, are stated to correspond to conceivable targets. The association of targets 
is improved and completed by corresponding calculations being carried out for 
Doppler velocities. Finally, the positions of the targets are calculated from a system 
of equations for the bistatically measured distances. 



